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ABSTRACT 

We reevaluate the hadronic part of the electromagnetic vacuum expectation 
value using the standard dispersion integral approach that utilizes the hadronic 
cross section measured in e~^e~ experiments as input. Previous analyses are based 
upon point-by-point trapezoidal integration which does not treat experimental er- 
rors in an optimal way. We use a technique that weights the experimental inputs by 
their stated uncertainties, includes correlations, and incorporates some refinements. 
We find the five-fiavor hadronic contribution to the fractional change in the electro- 
magnetic coupling constant at = m|, Aa(M|), to be 0.02752 ± 0.00046, which 
leads to a value of the electromagnetic coupling constant, a~^(M^) = 128.96±0.06. 
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1. Introduction 



At the current time, a large program of precise electroweak measurements is 
being conducted throughout the world. The object of this program is to test 
the electroweak Standard Model by comparing the measured values of a large set 
of electroweak observables with the predictions of the Minimal Standard Model 
(MSM). The Standard Model calculations have been performed to full one- loop 
accuracy and partial two-loop precision by a large community of researchers. In 
all of these calculations, it is necessary to evaluate the one-particle-irreducible 
contributions to the photon self-energy 11^^ (g^) or the related quantity U.!y^{q'^) = 
{'n.j^{q'^) — llj^{0))/q'^ at the Z mass scale — M|. These quantities are usually 
absorbed into the definition of the running electromagnetic coupling a{q'^), 

"^^'^ ^ i-[n;,(g2)°-n;,(o)]' 

where ao = 1/137.0359895(61) is the electromagnetic fine structure constant. This 
quantity is also represented as the fractional change in the electromagnetic coupling 
constant Aa, 

^«(^') - - - n;,(o). (2) 

Using analytic techniques and the optical theorem applied to the amplitude for 
s-channel Bhabha scattering, the quantity Aa has been related to the cross section 
for the process e+e" — > 7* — > all (atot) ^ follows [1], 



Aa(g^) = / ds^^^^RMs), (3) 

where Rtot{s) is the ratio of the total cross section to the (massless) muon pair cross 
section (t^^j_{s) — 47rQ;^(s)/3s at the center-of-mass energy ^/s. The cross section 
(Ttot is the physical cross section which has been corrected for initial state radiation. 
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The actual quantity measured in most experiments is discussed in Appendix A. 
It should be noted in passing that equation (3) is correct to all orders in ^-nd 
relies only upon the assumption that the real part of U^j is much larger than 
its imaginary part (the next-order correction is proportional to Im^rT-y-y/III-y-yp 
which is approximately 3xl0~^ at — M|). It is straightforward to evaluate 
equation (3) for the continuum leptonic cross sections [2]. In the limit that the 
scale is much larger than the square of the lepton mass mj, the contribution of 
the continuum leptonic cross sections is given by the following expression, 



.5 ^ ?1 

3 m« 



(4) 



The remaining contributions to Rtot consist of the continuum hadronic cross 
section and the = 1~ resonances and are labelled Rhad- Since the cross sec- 
tions for the resonances and low energy continuum are not accurately calculable 
from first principles, experimental inputs are used to evaluate their contributions 
equation (3). The contribution of open top quark production to the integral is accu- 
rately calculable and since the top quark mass is not known precisely, only the five 
flavor hadronic cross section is included in Rhad- The corresponding contribution 
to Aa{q'^) is therefore, 

a 7 2 

^(^hadiq^) = j ^^ g(g2^_g-) Rhad(s). (5) 

Equation (5) has been evaluated at the Z boson mass scale a number of times [3- 
8]. The most recent evaluations are by Martin and Zeppenfeld [6], Eidelman and 
Jegerlehner [7], and by Burkhardt and Pietrzyk [8] yield 

{ 0.02739 ± 0.00042, Reference 6 
0.0280 ± 0.0007, Reference 7 (6) 
0.0280 ± 0.0007, Reference 8. 

The authors of Reference 6 use perturbative QCD to parameterize the continuum 
Rhadi^) above y/s — 3 GeV and linear interpolation of measured data below that 
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point. The two-body final states 'k^ti~ and K^K^ are fit to parameterizations 
which include the p, c<j, and resonances. The remaining resonance contributions 
are calculated from an analytic expression which results from integrating a Breit- 
Wigner lineshape and depends upon the masses, widths, and leptonic widths of 
each resonance. The authors of Reference 7 use linear interpolation (trapezoidal 
integration) of measured data points to evaluate the continuum, 7r''"7r~, and K'^K~ 
contributions. Above ^/s — 40 GeV, they use perturbative QCD to evaluate Rhad- 
The contributions of the a;, 0, J/'^- family, and T-family resonances are included 
by integrating a Breit-Wigner lineshape. The authors of Reference 8 use smoothed 
averages of data to evaluate the continuum contribution, a parameterization to 
evaluate the tt+tt^ contribution, and the analytic expression to evaluate the con- 
tribution of the remaining resonances. 

This document reports on an evaluation of equation (5) which is performed 
in a somewhat different way from those listed above. In particular, the technique 
employed makes better use of the information provided by the various Rhad mea- 
surements, avoids some pitfalls inherent in the trapezoidal technique, and naturally 
provides an accurate estimate of the uncertainty on the result. We find 

Aahad{Ml) = 0.02752 ± 0.00046, 

which appears to be consistent with Refs. 6-8 within quoted errors. 

The result reported here updates an earlier result [9] which was more discrepant 
with Refs. 6-8. The updated value of AQ;/iarf(M|) is larger than the previous one 
by 8.6 X 10-4 for five reasons. The previous analysis used the six-flavor deflnition 
of Aa^ad which differs from the five-fiavor quantity by 0.6 x 10"^. A (hopefully) 
less controversial choice of ag(M^) shifts the result by —0.5 x 10^^. The fitting 
procedure used in the previous analysis was biased toward smaller Rhad values; 
correction of this problem gives a difference of 2.9 x 10"^. Small corrections to 
the analysis of the resonant contribution change the result by —0.1 x 10"^. But, 
the largest change is caused by the incorporation of a precise, new measurement of 
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Rhad near charm threshold which ahers the resuh by 5.8 x 10~^. Ahhough the net 
result is somewhat closer to those given above, a detailed comparison of the actual 
integrated cross section with one used in a trapezoidal integration (see Section 2.7) 
indicates that significant differences persist. 

2. The Analysis 

Any attempt to combine the results of many experiments is a perilous undertak- 
ing. Many different techniques and approaches have been used. Not all researchers 
have addressed all possible problems nor are systematic error estimates performed 
in uniform ways or to uniform standards. We therefore adopt some the techniques 
of the Particle Data Group [10]. Older measurements which are contradicted by 
newer, more precise work are excluded from the analysis. Parameter uncertainties 
that are extracted from fits with per degree of freedom (dof ) larger than one 
are rescaled by the factor \/~)^^Jdoj . 

2.1 Analysis Technique 

The experimental measurements of Rhadi^) ^-^e performed over limited regions 
of = y^. Typically, an experimental result consists of several points -R/j^d ~ 
Rhadi^i) measured at closely spaced energy points Wi. Each set of measurements 
is accompanied by a set of point-to-point uncertainties (statistical and systematic) 
(7j(ptp) and an overall normalization uncertainty (7(norm). Quite often, the point- 
to-point uncertainties are much smaller than the normalization uncertainty. A 
typical experimental result therefore consists of an accurately measured shape of less 
certain normalization. In this case, the values of the measured points are strongly 
intercorrelated. For future reference, we label these as Type I correlations. 

The normalization uncertainties usually incorporate purely detector-related 
effects, acceptance uncertainties, and uncertainties on radiative corrections and 
background corrections. The largest normalization uncertainties (15-20%) are as- 
sociated with the oldest measurements of Rhad iii the W — 1 — 5 GeV region. These 
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experiments typically had limited acceptance which when combined with a (com- 
mon) limited understanding of the event structure lead to large uncertainties in 
the overall detection efficiencies. The normalization errors associated with differ- 
ent sets of measurements performed at similar energies and times may be strongly 
correlated. These correlations are distinct from those discussed above (which must 
be present) and are labelled as Type II correlations. When combining the results 
of separate experiments, one must be careful to include the possible presence of 
Type II correlations in a conservative estimate of the overall experimental uncer- 
tainty. 

Most previous analyses of Aaf,ad evaluate various contributions to equation (5) 
by performing a trapezoidal integration with measured values of Rhad- Different 
data sets are combined by weighting nearby points by the quadrature sums of their 
point-to-point and normalization uncertainties (assuming that all points are uncor- 
related) . The effects of possible Type II correlations on the overall uncertainty are 
accounted for differently in different analyses. Eidelman and Jegerlehner [7] sum 
the uncertainties associated with each point linearly. Burkhardt and Pietrzyk [8] 
and most of the earlier analyses assign typical normalization uncertainties to vari- 
ous intervals in W and sum the corresponding uncertainties on Aa^ad in quadra- 
ture. The use of trapezoidal integration has two advantages: it is unbiased by 
human prejudice about the functional form of Rhad{s), and it would automatically 
account for undiscovered resonances which are broad as compared with the spacing 
of measurements. Unfortunately, this technique also has a serious shortcoming: it 
ignores the Type I correlations present in each data set. 

Treating the combined (normalization and point-to-point) uncertainties on the 
points in each set as independent loses the (often precise) shape information asso- 
ciated with the set. Two examples of the loss of shape information are illustrated 
Figure 1. In part a), a data set with small point-to-point errors (shown as solid 
dots) and a large normalization uncertainty (illustrated to the right of the data) 
is combined with a single precise measurement (shown as the open dot) . The sta- 
tistical averaging procedure used in the trapezoidal integrations would yield the 
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function shown as the sohd curve. The shape defined by the sohd dots would be 
distorted near the single precise point and the accurate normalization information 
contained in the single measurement would be ignored. An more optimal procedure 
would use the shape information provided by the solid dots and the normalization 
information provided by the open dot yielding the dashed curve. 
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Method \ 
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Estimate 
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Figure 1. Two examples of the shape infor- 
mation loss inherent in the averaging procedures 
used by trapezoidal analyses. 

Part b) of Figure 1 shows the result of combining two partly overlapping sets 
which have small point-to-point uncertainties and large normalization uncertainties 
(shown as open and solid dots, respectively). In the region of overlap, the sets define 
a consistent shape but differ in normalization. An optimal averaging procedure 
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would average the normalizations and produce the dashed curve. The procedure 
adopted as part of the trapezoidal analyses would yield the solid curve which agrees 
with the dashed one only in the region of overlap and does not preserve the shape 
determined by the data sets. The trapezoidal analyses described in References 4, 5, 
and 7 are checked by first integrating individual data sets and then by averaging the 
integrals. While it might appear that this procedure preserves shape information, 
the actual averaging of the integrals can be carried out only in energy intervals 
where the data sets overlap. The net result therefore looks much like the solid curve 
in part b). It is not surprising that consistent results were obtained. Optimal use of 
the shape information can occur only in techniques that allow the normalizations 
of the data sets to vary. The consequences of these examples will become clearer 
in Section 2.7. 

We incorporate correlations into the analysis by fitting the data to an appro- 
priate functional form Rfit{s; a^) where are the parameters of the function. In 
the absence of undiscovered resonances, Rhad can be described by a continuous 
function. A fit has the virtue that measurements can be weighted by their 
experimental errors and correlations are straightforward to include. The previous 
version of this analysis used a non-diagonal definition of constructed from the 
covariance matrix Eij = (^ARj^^^ARj^^^^ of the measured points R\^^- Unfor- 
tunately, it has been shown that if the off-diagonal elements Eij scale with the 
measured values of -R^^^ or Rj^^^, the resulting fit will be biased to smaller values 
of Rhad [11]- The bias that resulted to our previous analysis [9] from the appfica- 
tion of the incorrect technique was approximately 39% of the uncertainty on the 
final result. We avoid the bias by defining as follows. 



where R^g^^ is the value of Rhad measured at energy Sj, ai = ai{noTm) / Rj^^ is the 
fractional normalization uncertainty associated with the i^^ measurement, and Xj 
are fit parameters which are constrained to have zero mean and unit width. This 



[^M - (1 + Aja,:)-R/a(8i; at)] 
<T?(ptp) 



j 



(7) 
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form preserves shape information and propagates the normaUzation uncertainties 
into the parameters of the function Rfn- For each fit, two choices of the parame- 
ters A,- are investigated. In the first separate normahzation parameter A, 
is assigned to each data set. This choice incorporates Type I correlations only and 
makes no assumptions about correlations between experiments. In the second case, 
the normalizations of experiments of similar age and energy region are assumed to 
be 100% correlated. A separate normahzation parameter is assigned to each cor- 
related group instead of each set of measurements. This choice includes the effects 
of Type I and Type II correlations, produces larger error estimates (a consequence 
of including the Type II correlations), and is the one quoted as the official result. 
The difference in Aahad resulting from the two weighting schemes is included in 
the parameterization uncertainty discussed below. 

Equation (5) is evaluated by performing a Simpson's Rule integration using 
the function Rfn and the best estimate of the parameters. The parameter un- 
certainties 5ak reflect the point-to-point and normalization uncertainties to some 
extent. Unfortunately, the process of fitting a large number of measurements with 
a function of a smaller number of parameters necessarily involves some loss of in- 
formation. The resulting uncertainty on the fitting function at some point W is 
usually smaller than the uncertainties on nearby data points. If we add a priori 
information to the problem by choosing a physically motivated fitting function, 
the information contained in the parameter error matrix may be appropriate. To 
understand this problem better, we evaluate the uncertainty on Aa^adi^z) 
two techniques. In the first, the parameter uncertainties are propagated to the 
calculated value of Aahad{M^) using the following expression which is valid for 
any function of the parameters, 

A2rA^ ^ djAahad) ^ djAahad) 

d [Aahad)exp = V — ^ Eki , (8) 

k,l 

where the derivatives are calculated numerically and E^i — {Sa^Sai) is the param- 
eter error matrix that is extracted from the fitting procedure. The second error 
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estimate is performed by constructing a large ensemble of data sets by shifting the 
measured data points i?J^^^(meas) as follows, 

i?Ld(set j) = Riadi^eas) + ffj^'^^iivtv) + (9) 

where the factors /y are Gaussian-distributed random numbers of unit variance. 
The entire fitting and integration procedure is then applied to each member of the 
ensemble. The uncertainty on Aafiad{M^) is determined from the central 68.3% 
of the ensemble distribution. 

The use of a fitting function has the problem that one may introduce bias 
through the choice of parameterization. We attempt to evaluate this effect by 
varying the parameterizations as much as ingenuity and computer time allow. The 
quoted contributions to Aahad{M^) are those corresponding to the best fits. Each 
contribution is assigned a parameterization uncertainty S{Aa}iad)param based upon 
the spread of results corresponding to reasonable fits. The parameterization un- 
certainty also includes a contribution from the difference observed in the two 
weighting schemes. 

2.2 The Data 

The approach to the evaluation of equation (5) is driven by the form of the 
data themselves. The total hadronic cross section can be decomposed into four 
pieces: the hadronic continuum above W = y/s = 1 GeV, the charged two-body 
final states tt+tt^ and K^K~~ from their respective thresholds to 2.6 GeV, and 
hadronic resonances (excluding charged two-body final states). Since equation (5) 
is linear in the hadronic cross section, we decompose Aajiad as follows, 

Aauadiq^) = A<^^(g2) + Aa^f {q^) + AaL?^?') + ^^Mdi^'l W 

where the four terms on the right-hand side correspond to the four pieces of the 
hadronic cross section. 
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The rationale for this decomposition is as follows. The region below W = 
1 GeV is dominated by the p, u, and resonances. The electromagnetic form 
factors for the processes e~^e~ — > 7i^7i~ [12-19] and e+e' [19-23] are 

measured well from threshold to ~ 2 GeV. Resonances do not account for all of 
the tt+tt" and K'^K~ cross section in this region. On the other hand, essentially 
all other two-body and three-body final states are associated with the resonances. 
Measurements of three-pion final states near W — 1 GeV [24] show the non- 
resonant portion to be consistent with zero. Similarly, measurements of various 
two-body final states such as Kj^Kg show small non-resonant cross sections [21]. 
The cross sections for four-pion final states become significant above 1 GeV but 
are small below that energy [25]. The 772 experiment [26] at the ADONE storage 
ring at Frascati has measured the hadronic cross section ratio for three or more 
hadron final states, R^^^ from W = 1.42 GeV to W ^ 3.09 GeV. They have also 
presented several points from 1 GeV to 1.4 GeV that are composed of various 
multipion cross sections from Novosibirsk and Orsay [25,24,27] and are claimed to 
approximate R^^^- Measurements beginning at 1^ = 2.6 GeV by the MARK I [28], 
DASP [29], PLUTO [30], and Crystal Ball [31] Collaborations claim to measure the 
entire cross section. We therefore conclude that Rhad is well approximated below 
Wi = 2.6 GeV by a sum of the 7r"'"7r~ and K'^K~ contributions from threshold to 
Wi (where they are much smaller than Rj^^j)', the Rj^^^ measurements from 1 GeV 
to Wi; and the p, lu and (p resonances where the hadronic widths are adjusted 
to remove the 7r''"7r~ and K^K~ final states that are already included explicitly. 
Note that the several broad e+e~ resonances between the 0(1020) and W — 2 GeV 
are implicitly contained in the two-body or R^^^ categories. Since the 7r"'"7r~ and 
K'^K~ cross sections are very small at Wi, the R^^^ and total continuum Rhad 
measurements should be continuous at this point. 

At center-of-mass energies larger than Wi, many measurements of the hadronic 
continuum and resonances exist. The only precise measurement in the region 
between 2.6 GeV and 5.0 GeV is a single data point just below charm threshold at 
W — 3.670 GeV by the Crystal Ball Collaboration [31]. This measurement has a 
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normalization uncertainty of 7%. Since the next most precise measurements in the 
region below 5 GeV have normalization uncertainties of 15%, this measurement 
represents an important constraint on the magnitude of the cross section in the 
entire region. The region above charm threshold from W — 3.77 GeV to = 
5.0 GeV is comphcated and not well measured. The MARK I, DASP, PLUTO 
and Crystal Ball Collaborations all observe an enhancement beyond the expected 
threshold shape. The DASP data show three resolved resonances. The MARK I 
and PLUTO data are consistent with the DASP data but do not cleanly resolve the 
resonances. The Crystal Ball measurements are somewhat smaller than the older 
ones and do not resolve the second resonance (which appears as a broad shoulder). 
We choose to follow the Particle Data Group and recognize the DASP resonances: 
-0(4040), '0(4160), and -0(4415). The ip family therefore consists of six states. 

Between 5 GeV and 10.4 GeV, the MARK I, DASP, PLUTO, Crystal Ball [32], 
LENA [33], CLEO [34], CUSB [35], and DESY-Heidelberg [36] Collaborations have 
published Rhad measurements which are plotted in Figure 2. The error bars in- 
clude only point-to-point uncertainties. The recently published Crystal Ball mea- 
surements have a systematic normalization uncertainty of 5.2%. The other mea- 
surements have normalization uncertainties in the range 6.8-10%. The data are 
also compared with the recent QCD prediction of Chetyrkin and Kuhn [37] which 
includes quark mass effects. At 1^ = 5 GeV, the MARK I data are consistent with 
other measurements. As W increases, they show a systematic increase in Rhad 
and suggest the presence of a structure near 6.6 GeV. Including the quoted 10% 
normalization uncertainty, the MARK I data are larger than the more precise mea- 
surements by approximately two standard deviations. The reader is reminded that 
first generation detectors like MARK I, DASP, and PLUTO were small acceptance 
devices that necessarily involved large acceptance corrections without the benefit of 
good event structure modelling. After acceptance corrections and a r-lepton sub- 
traction, the MARK I group observed that two-charged-prong events constituted 
nearly 20% of the hadronic cross section of i? at 1^ = 7 GeV. This is about 1.5 
times the two-prong rate due to t'^t~ final states and three times the rate that is 
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predicted [38] by the JETSET 7.3 Monte Carlo program [39]. While this may not 
be wrong, we choose to exclude data from the first generation experiments when 
more modern results are available. Such data are available above charm threshold. 
Unfortunately, we are constrained to use very old continuum measurements below 
charm threshold. 
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Figure 2. The Rhad measurements of the MARK I [28] (sohd 
squares), PLUTO [30] (open triangles). Crystal Ball [32] (open di- 
amonds), LENA [33] (open squares), DASP [29] (open inverted 
triangle), and DESY- Heidelberg [36] (open circle) Collaborations 
in the region between W = 5 GeV and W = 9.4 GeV. The error 
bars include point-to-point uncertainties only. A recent QCD cal- 
culation [37] which includes quark mass effects is shown as a solid 
line for Qs(M|) = 0.125. 

The Particle Data Group lists six T family resonances between 9.4 GeV and 
11 GeV. All are included in the resonance contribution. 

Above b-quark threshold, a number of Rhad measurements have been carried 
out by the PEP and PETRA experiments [40-45]. However at energies above 
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W = 3A GcV, Z-7 interference becomes significant. We therefore use only those 
measurements in the region W < 34 GeV where the correction for electroweak 
interference is less than 1%. 

We expect that Rhad is well described by perturbative QCD in the region 
above b-quark threshold. This implies that the world average value of the strong 
coupling constant a;s(M|) compiled by the Particle Data Group [10] provides a 
precise measurement of Rhad aX W = Mz- Since possible anomalies in the Z 
lineshape would bias the determination of as{M'^) from the lineshape parameters, 
we exclude the Z lineshape information from the Particle Data Group average. 
Additionally, since we explicitly include the PEP/PETRA Rhad measurements in 
our fit (which uses perturbative QCD to describe the PEP/PETRA energy region), 
they are also excluded from the PDG average yielding the following value. 



q;s(M|) = 0.116 ±0.005. 



(11) 



To convert as(M|) into a determination oi Rhadi^z), we use the third-order QCD 
expression [46], 



RQCD{s)^3j2Q}Pf 



(3 - 



2 



/ 



(12) 




+ ri 




TT 



+ r2 




where: Qf is the final state fermion charge, (3f = ^1 — imj/s is the fermion veloc- 
ity in the e~^e~ center-of-mass frame (my is the fermion mass), and the coefficients 
are functions of the number of active flavors Nf, 



n = 1.9857 - 0.1153A^/ 



r2 = - 6.6368 - 1.2002A^/ - 0.0052A^| - 1.2395 



(13) 
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The resulting value of R^adi^z) is, 

i?/iad(Mz) = 3.807 ±0.006. (14) 

The following three sections of this chapter describe the evaluation of: the 
continuum contribution Aq;^°^*^, the contributions of the charged two-body final 
states Aa^^J" and Aa^^^ , and the resonance contribution /S.a^^^. 

2.3 The Hadronic Continuum 

The first step in the evaluation of equation (5) for the hadronic continuum 
is to formulate a suitable (piecewise-continuous) parameterization Rfn{s]a]^). We 
choose to use the perturbative QCD expression given in equation (12) in the region 

W > lb GeV and an empirical parameterization in the region 1 GeV< W < 
15 GeV. In the high energy region, the only free parameter is as{M'^) which is 
evolved to other scales numerically using the Runge-Kutta method applied to the 
order-ctf renormalization group equation [47]. 

In the portions of the low energy region that are measured well, polynomials are 
used to parameterize RhadiW). To ensure that the function is continuous across 
several points Wa-, the polynomials are constructed in Xa = W — Wa and the zeroth 
order terms are excluded, 

n 

P«(x)^J]<4, (15) 
1=1 

where a is a label to distinguish different regions. Separate polynomials are used 
to describe the following regions: 1 GeV< W < 1.9 GcV (labelled region s), 
1.9 GeV< W < 3.6 GeV (labelled region c), and 5.0 GeV< W < 10.4 GeV (labelled 
region b). Although a single, large-order polynomial is adequate to describe the 
data between W = 1 GeV and charm threshold at 3.6 GeV, the data show a 
distinct shape change near 1^ = 1.9 GeV (where the four-pion cross section is 
becoming small) . It was possible to obtain better fits by introducing an additional 



15 



polynomial to describe the region from 1 GeV to 1.9 GeV. A comparison of the two 
possible forms is used to assess the parameterization sensitivity of the final result. 

Since there are no measurements of the continuum Rhad the b-quark and c- 
quark threshold regions (published measurements include a mixture of continuum 
and resonances), it is necessary to extrapolate the form of Rhad from 3.6 GeV to 
5.0 GeV and from 10.4 GeV to 15 GeV with functions that are physically motivated. 
In the case of the charm threshold region, the DASP Collaboration has pubhshed 
(in graphical form) the shape of the continuum that was preferred by their fit to the 
■?/'(4040), ■?/'(4160), and ■(/^(4415) resonances. The function which characterizes the 
shape of the threshold, fDASPi^)^ does not increase as sharply as the free-quark 
threshold factor /?(3— /?^)/2 but increases more rapidly than the threshold factor 
for pointlike scalar particles. To construct the function Rfn, all three possibilities 
are used for the c-quark threshold and the two extreme possibilities are used for 
the b-quark threshold. 



where the c- and b-quark masses are taken to be the D and B meson masses, 
respectively. The actual size of the charm-associated step in Rhad: is left as 
a free parameter. The size of the bottom-associated step in Rhad is constrained to 

be the difference between the value of the fit function at 1^ = 10.4 GeV and the 



value of the QCD portion at W = 15 GeV, AR,, = Rqcd{^^) - RfuilOA). 
The actual form of the fitting function is given by the following expression. 




(16) 



' Ro + P^^{W -1.0), 1<W^<1.9 

R fitil. 9) + P'ji^^iW -1.9), 1.9<VF<3.6 

Rfit{S.6) + ARJ.iW), 3.6 < 1^ < 5.0 

Rfit{5.0) + PI^^{W -5.0), 5.0 <W^< 10.4 

RfitilO.4) + ARhhiW), 10.4 <W< 15.0 

. Rqcd{W), 15 <W 



(17) 
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where Rq, the value of Rhad at = 1 GcV, is a free parameter and the or- 
der of the polynomials is varied from 1 to 7. The is constructed from equa- 
tion (7) assuming that normalization uncertainties are completely correlated in 
four groups: the 20% uncertainties of the lowest energy measurements [26-27] 
(1.0 GeV< W < 3.09 GeV), the 15-20% uncertainties of the MARK I, DASP, and 
PLUTO measurements [28-30] (2.6 GeV< W < 4.9 GeV), the 5-10% uncertainties 
of the measurements [32-36] between charm and bottom thresholds (the Crystal 
Ball measurement at 3.670 GeV is treated as a member of the higher-energy Crystal 
Ball set), and the 1.7-7.0% uncertainties of the PEP and PETRA experiments [40- 
45] above bottom threshold. Each fit is repeated with a separate normalization 
parameter for the 20 sets of data in the analysis. 

The data are corrected for electroweak interference and incomplete vacuum 
polarization corrections (see Appendix A) before the fitting procedure is applied. In 
the course of varying the orders of the polynomials and the number of normalization 
parameters, the number of free parameters varies from 12 to 44. The fit quality 
does not improve substantially when the number of parameters exceeds 14. The 
data and the result of the fit used to to evaluate the central value of AaJ^^^ are 
shown in Fig. 3. The error bars include the point-to-point and the normalization 
uncertainties. The fit quality is reasonable (x^/dof = 110/100). 
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Figure 3. The continuum Rhad measurements including normalization uncer- 
tainties. The entries in the region below charm threshold consist of a compilation 
of low energy exclusive cross sections [25,24,27] (solid inverted tranglcs) and the 
measurements of the 772 [26] (solid dots), Mark I [28] (open diamonds), DASP [29] 
(X's), Crystal Bah [31] (solid square), and PLUTO [30] (soHd diamond) Collabo- 
rations. The entries in the region between charm and bottom thresholds are the 
measurements of Crystal Ball [32] (open diamonds), PLUTO [30] (open triangles), 
LENA [33] (sohd squares), DASP [29] (diamond-X overlay), DESY-Heidelberg [36] 
(square-X overlay), CUSB [35] (solid dot), and CLEO [34] (sohd inverted triangle) 
Collaborations. The entries in the region between above bottom threshold and be- 
low the Z pole are the measurements of CELLO [42] (open diamonds), PLUTO [30] 
(open triangles), JADE [33] (open squares), Mark- J [44] (open inverted triangles), 
TASSO [45] (circlc-X overlay), HRS [40] (open circle), and MAC [41] (X) Collabo- 
rations. The fit used to to evaluate the central value of Aq;^°'^' is shown as the solid 
curve. 
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The various hypotheses for Rfn are used to evaluate the integral in equation (5) 
from So = 1 GeV^ to oo = 10^ GeV^. Although the singularity in the integrand is 
formally well controlled, digital computers are famous for their inability to under- 
stand formalities. We have therefore recast equation (5) into a form which is more 
suitable for electronic evaluation, 



dRfit 



0\ 



ds 



In 



- so' 

So . 

+ A 



So 

oo 



q 



2- A 



s{s — q^) 
s{s — (p) 



(18) 



where we have assumed that Rfn is well approximated by a linear expansion over 
the interval — IS. < s < -\- iS. (ia practice, we use A = 0.5 GeV^). The 
evaluation of equation (18) requires that as be evolved to scales larger than the 
t-quark mass. For this purpose, the top quark mass is assumed to be 172.3 GeV 
which is the MS mass corresponding to a pole mass of 180 GeV. 

The contribution of the hadronic continuum to Aa/ja^^ (M^) is found to be fairly 
insensitive to the form of Rfn and the number of normalization parameters used. 
The central value of AQ;/jod(-^f ) corresponds to the best estimate of the parameters 
of the function which uses: the DASP shape for the c-quark-threshold, the free- 
quark shape for the b-quark-threshold, the values (2,3,3) for {Nf„Nc,Ns) and four 
normalization parameters. The maximum deviation from this value occurs when 
Ni, = 1 and four (instead of 20) normalization parameters are used (the deviation is 
insensitive to the choice of threshold functions). The size of the maximum deviation 
is taken as an estimate of the parameterization uncertainty. The experimental 
uncertainty given by equation (8) is found to be in excellent agreement with the 
estimate derived from an ensemble of 500 fluctuated data sets. The resulting 
contribution to Aafiad{M'^) is 

Aal°^^{Ml) = 0.022106 ± 0.000366(exp) ± 0.000196(param). (19) 
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Figure 4. The uncertainty on the integrand of the dis- 
persion integral (integrated over W rather than s) in arbi- 
trary units. The dashed curve shows the uncertainty before 
the Crystal Ball data point is included in the fit and the 
solid curve shows the uncertainty after its inclusion. 

This result differs from our previous result [9] by +0.000678. Most of the differ- 
ence is caused by inclusion of Crystal Ball data point at 3.670 GeV (+0.000575). 
The remaining difference is due to the use of the five-flavor definition of A.afiad 
(+0.000059), a change in the value of q;s(M|) used as input (-0.000051), and 
the change to the unbiased fitting technique (+0.000095). The inclusion of the 
Crystal Ball point pulls the fit to somewhat larger values of Rhad and substan- 
tially constrains the normalization in the charm threshold region. The Mark 11 
and 772 data span a large energy region and constrain the shape of Rfit{W) down 
to W — lA GeV. The effect of the single precise point is therefore propagated 
to to smaller energies. This type of effect is illustrated in Figure 1(a) and is 
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demonstrated in Figure 4 which displays the uncertainty on the integrand of the 
W^-space dispersion integral in arbitrary units [48]. The uncertainty is calculated 
using equation (8) (with Aahad replaced by Rjit) to estimate the uncertainty on 
RfitiW) at each energy point. The dashed curve shows the uncertainty before the 
Crystal Ball data point is included in the fit and the solid curve shows the uncer- 
tainty after its inclusion. Note that the overall uncertainty on Aq;^°^*^ is dominated 
by the poor precision of the data in the 1 GeV to 3.5 GeV region. 

2.4 The tt+tt" and K^K~ Final States 

The processes e+e^ tt+tt^ and e^e~ K^K^ are described by the elec- 
tromagnetic form factors, -^^(s) and Fk{s), which are related to the hadronic cross 
section ratio Rhad ^ot each process as follows, 

where (3^^ and (3k are the velocities of the final state particles in the e'^e~ center- 
of-mass frame. It is clear that measurements of the form factors are equivalent to 
measurements of Rhad- 

Measurements of the square of the pion form factor l-FVp have been performed 
by the OLYA [12], CMD [12], TOF [14], NA7 [13], /xtt [17], MEA [19], M2N [15], 
DM1 [16], and DM2 [18] Collaborations and are shown in Fig. 5. The error bars 
include the normalization uncertainties which range from about 2% in the region 
around the (dominant) p resonance to about 12% at ~ 2 GeV. 

The data are first corrected for incomplete vacuum polarization corrections as 
described in Appendix A. They are then fit to a function which is a sum of the 
Gounaris-Sakurai form [49] used by Kinoshita, Nizic, and Okamoto [50] and three 
resonances, 

p (A^ Ai-mlA2 A Bne^^-ml 

^""^'^ A, + An^ + /(,) + A. , - m2 + ^rnnVn ' ^ ^ 

where: Ai and A2 are free parameters; ttItt is the pion mass; q and f{s) are defined 
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as follows, 

q = ^ s/4 - ml 
f{s) = - - - + -^In 

and where m^, Tn, -Bn, and Cn are the mass, width, amplitude, and phase of 
each resonance. The mass and width of the first resonance are set to those of the 
a;(782). All other parameters (12 in total) are allowed to vary. The function 
is constructed assuming that all normalization uncertainties are 100% correlated 
(one normalization parameter) and that the normalizations are uncorrelated (seven 
normalization parameters). As in the case of the continuum, the two fits give 
nearly identical results but the error estimate is larger when only one normalization 
parameter is used. The result of the single-normalization-parameter-fit is shown 
as a solid hne in Fig. 5. The fit preferred a resonance of width 0.44 GeV at mass 
1.15 GeV and a second resonance of width 0.18 GeV at mass 1.71 GeV. The fit 
quality is found to be good (x^/dof = 138.3/127). 

To evaluate the sensitivity of the result to the parameterization, the complete 
function used by the authors of Ref. 50 was also fit to the data. This function did 
not fit the newest (large W) data from DM2 as well as our chosen form (x^/dof = 
201.5/132). Both functions were used to evaluate equation (5) from s = Am^ 
to s = 4 GeV^ (where I-FttP is measured to be very small). We find the 7r"'"7r~ 
contribution to ^a^adi^'z) 

AaJ^^j" (M|) = 0.003240 ± 0.000057(exp) ± 0.000169(param). (23) 

The two techniques for the estimation of the experimental uncertainty (discussed 
in Section 2.1) yield consistent results. 

The result given in equation (23) differs from our previous result [9] by 
-1-0.000153. The difference is due entirely to the use of the unbiased fitting tech- 
nique and represents the largest problem found with the older technique. 
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Figure 5. Measurements of \F^{W)\^ by the OLYA [12] (solid dots), CMD [12] 
(open diamonds), TOF [14] (solid triangles), NA7 [13] (open squares), /ivr [17] (solid 
squares), MEA [19] (solid diamonds), DM1 [16] (open triangles), and DM2 [18] (open 
circles) Collaborations are compared with the best fit which is shown as a solid line. 
The error bars include normalization uncertainties. 

Measurements of the square of the kaon form factor \Fj{\'^ have been performed 
by the OLYA [20], CMD [21], MEA [19], DM1 [22], and DM2 [23] Collaborations 
and are shown in Fig. 6. The data span the 0(1020) resonance and continue to 
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W = 1.8 GcV where -R^^^^ is less than 0.01. The normahzation uncertainty 
on the CMD measurements is 6%. The other groups do not report normalization 
uncertainties. Early IF^]"^ measurements suffered from the same problem of un- 
reported normalization uncertainties. A bit of historical research shows that the 
normalization uncertainties were usually not included in the point-to-point errors. 
We therefore arbitrarily assign a 20% systematic normalization uncertainty to all 
unreported cases. The data and total uncertainties are shown in Fig. 6. 




5-95 W (GeV) 7854A4 



Figure 6. Measurements of \Fk(W)\^ by the OLYA [20] (solid dots), CMD [21] 
(open diamonds), MEA [19] (open squares), DM1 [22] (open triangles), and DM2 [23] 
(open circles) Collaborations are compared with the best fit which is shown as a solid 
line. The error bars include normalization uncertainties. 

The data are fit to a function which is a sum of a Breit-Wigner resonance with 
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an energy- dependent width for the (j) and four resonances, 

A-, ^ R 

where: Ai is the amphtude of the 0; is the mass of the 0(1020); mm r„, Bm 
and Cyj are the mass, width, amphtude, and phase of the resonances. The energy- 
dependent width r^(s) is assumed to consist of contributions from the K'^K~ , 
K^Kg, and Stt final states. 



I3l(s) (3^(s) 
0.497 ,,+ ^ +0.347-^°^ ^ 



+ 0.156Gi{s)\, (25) 



where: is the nominal value [10] of the (j) width, (3+{s) — 4m|.+ /s is the 

velocity of the charged kaon, Po{s) — -y/l — Am^o/s is the velocity of the neutral 
kaon, and ^3^(5) is a function which is normalized to unity at s = and is 
proportional to the decay rate for (p ^ Sir assuming pir dominance [51]. 

The masses and widths of the first two resonances were set to those of the p(770) 
and u;(782). Following the procedure of Ref. 23, the amplitude ratios Bi/Ai and 
B2/A1 were constrained to the measured values and the phases were set to zero. 
The mass, width, and amplitude of the were allowed to vary. The masses, widths, 
amplitudes, and phases of two larger mass resonances were free parameters. The 

function was constructed with the assumptions that all normalization uncer- 
tainties are 100% correlated (one normalization parameter) and the normalization 
uncertainties are uncorrelated (five normalization parameters). The fit was 

the only instance for which the different assumptions about the correlation of the 
normalizations yielded noticeably different fit results. In this case, the assumption 
that the normalizations are uncorrelated (five normalization parameters) produced 
a substantially better fit to the data (x^/dof = 48.9/44) than did the assumption 
that they are correlated (x^/dof = 73.6/48). The better fit is plotted as a sohd 
line in Fig. 6. The fit preferred a resonance of width 0.17 GeV at mass 1.35 GeV 
and a second resonance of width 0.24 GeV at mass 1.68 GeV. 
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To evaluate the sensitivity of the result to the parameterization, a second fit 
was performed with the amplitudes and phases of the p and oo allowed to vary as 
free parameters. No appreciable differences differences from the first pair of fits 
were observed. Evaluating equation (5) from s — 4m|.+ to s = 3.24 GeV^, we find 
the K^K~ contribution to Act/jadl-^I) ^® 

^(^hlf i^z) ^ 0.000356 ± 0.000032(exp) ± 0. 000030 (param) (26) 

where the parameterization uncertainty refiects the difference obtained from the 
two definitions. The two techniques for the estimation of the experimental 
uncertainty (discussed in Section 2.1) yield consistent results in this case. 

2.5 The Resonances 

The resonances comprise the remaining portion of the total e'^e~ cross section. 

The total cross section for each resonance can be represented by a relativistic Breit- 
Wigner form with an energy-dependent total width [52], 

_ 12. y-sr^erfsis) 

where: m, Fge, and Ttot are the mass, electronic width, and energy-dependent total 
width of the resonance; and F fg is the energy-dependent width corresponding to the 
final states considered in the analysis. Note that the electronic widths are physical 
widths (not corrected for vacuum polarization effects) . In order to incorporate the 
Breit-Wigner cross section described by equation (27) into equation (5), it must 
be scaled to the electromagnetic point cross section, cfhh{s) — 47rQ;^(s)/3s, yielding 
the following expression, 

2x "Or^ [ , (tresis) 



= J (28) 



4m: 



2 



which has the slightly unpleasant feature that it incorporates a{s), the quantity 
that we are attempting to evaluate, into the integrand. To avoid this problem, we 
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use the Aahadis) parameterization given in Ref. 4 to generate a first-order estimate 
of a{s) for use in equation (28). Note that equation (28) is often written with a{s) 
replaced by ao. This is correct only if the cross section aj-es is replaced by the 
tree-level one, a^^g — ares • a.Q/a,'^{s). The factor aQ/a^{s) is often absorbed into 
equation (27) by defining the tree-level electronic width Fgg = Fee • aQ/a'^{m'^). 

Equation (28) is evaluated for the a;(782), 0(1020), V'-family, and T-family 
resonances by performing a Simpson's rule integration over the interval m — QOTtot 
to m + QOTtot (the lower limit of the uj integration is the threshold for Sir decay). 
The energy-dependent total widths of the ip and T resonances are assumed to scale 
as y/s, 

Ttotis) = ^Ttotim), (29) 
m 

where m is the mass of the resonance and rtoi(m) is the nominal value of the 
width. All ip and T final states are included in the resonance contribution [rjs(s) = 
Ttot{s)]- The energy-dependent total width of the 0(1020) is given by equation (25). 
The width rjs(s) for the is adjusted to exclude the K^K~ final state (discussed 
in Section 2.4). The energy-dependent total width of the c<j(782) is given by the 
following expression which assumes that all final states are tt+tt^, tt'^'J, or 7r+7r~7r°, 



m. 



0.022#4 + 0.085- (^-"^'/^^^ 



P^i^l) ■ (1 - m2/m2 )3^ 



+ Q.893Gt{s)}, (30) 



where: rriui is the mass of the u, is the nominal value [10] of the u width, 
P-k(s) — a/I — 4m2 / s is the velocity of the charged pion, and G'^^{s) is a function 
which is normalized to unity at s = and is proportional to the decay rate for 
a; — > Stt assuming a constant matrix element (phase space weighting). The width 
rjs(s) for the cu is adjusted to exclude the 7r"'"7r~ final states which are included in 
the l-FyrP contribution. 

The masses and widths used to evaluate equation (28) are taken from the 1994 
Review of Particle Properties [10]. The Particle Data Group does not apply a con- 
sistent set of definitions to the parameters of all resonances. The electronic widths 
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of the ip and T families are defined to be the physical ones and are derived from 
fits performed by the PDG itself. The electronic widths of the uj and resonances 
are determined from measurements of the total widths and electronic branching 
fractions Bgg. In both cases, the total widths are the correct physical ones. The 
average value of See((^) is dominated by peak cross section measurements of the 
CMD [53] and ND [54] Collaborations which are corrected (partly) for vacuum 
polarization effects and lead to a determination of r^^{u;). The case of the is less 
clear. Of the three most precise measurements of Bee{4>), those of the DM1 [55] 
and OLYA [56] Collaborations are not corrected for vacuum polarization effects 
and lead to a determination of ree((;i). The most precise measurement is a later 
OLYA result which has about the same precision as the combination of the two 
preceeding results but is reported in an unpublished preprint which is no longer 
available for inspection. The result may (or may not) be corrected for vacuum 
polarization effects. We make the assumption that the RPP value of ree(0) is 
the physical one. This assumption cannot be wrong by more than one half of the 
total vacuum polarization correction (1.6%) which we include in the uncertainty 
on ree(0). 

The leptonic widths are corrected for incomplete vacuum polarization correc- 
tion to the normalizing cross sections (see Appendix A) before equation (28) is 
evaluated. The results are listed in Table 1 along with those derived in Sections 2.3 
and 2.4. The experimental uncertainties are evaluated by assuming that the un- 
certainties on the masses, total widths, electronic widths, and relevant branching 
ratios are uncorrelated. The parameterization uncertainties are evaluated by re- 
peating the calculation with a constant-width, constant-mass Breit-Wigner cross 
section. 
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2.6 Final Result 

The various contributions to Aa^adiM'^) are summarized and summed in Ta- 
ble 1. The resulting value is 

Aahad{Ml) = 0.02752 ± 0.00046. (31) 

Including the leptonic contribution, we find q;-1(M|) to be, 

q;-1(m|) = 128.96 ± 0.06, (32) 

where the uncertainties on the lepton masses contribute negligibly to the total 
uncertainty. This result differs by one of its standard deviations from the (common) 
result given in References 7 and 8 and it differs by 0.3 standard deviations from the 
result given in Reference 6. However, since the different analyses make use of many 
of the same inputs, the results are not independent measurements of Aafiad{M'^) 
but reflect differences in assumptions and technique. 



Table 1: Summary of the various contributions to Aahad- 



Contribution 


W Region (GeV) 


AahadiMl) 


^{^(^had)exp 


^{^(^had)param 


Continuum 


l.O-oo 


0.022106 


0.000366 


0.000196 




0.280-2.0 


0.003240 


0.000057 


0.000169 


K+R- 


0.987-1.8 


0.000356 


0.000032 


0.000030 


Resonances 




0.000307 


0.000010 


0.000003 


II 




0.000296 


0.000012 


0.000004 


II 


ip (6 states) 


0.001101 


0.000059 


0.000023 


II 


T (6 states) 


0.000118 


0.000005 


0.000003 


Total 




0.02752 


0.00038 


0.00026 



Doesn't include tt+tt final states. 
Doesn't include K'^K~ final states. 
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2.7 Detailed Comparison With Reference 7 



The result of Eidelman and Jegerlehner [7] (henceforth E-J) is based almost 
entirely upon the trapezoidal integration of locally averaged data points. Only the 
narrow resonances are treated parametrically. E-J have published their composite 
compilation of the function Rhad{^) in a series of figures and include a detailed 
breakdown of the contributions of various energy intervals to AQ;/jad(-^f )■ Since 
the E-J compilation excludes narrow resonances, we construct the function Rsum 
to include the same final states, 



RsumiW) = RfuiW) + ^\F^{W)\^Pl + ^\Fk{W)\^PI + ^ al,,{W), (33) 

1=1 

where the sum includes the a;(782), 0(1020), ^(4040), ^(4160), and ^(4415) res- 
onances. A comparison of their Rhad compilation {Rfi^^) with Rsum in the region 
= 1 — 50 GeV is shown in Figures 7 and 8. The R^^^ compilation is shown as 
the solid curve in both figures. The dashed curve in Figure 7 shows Rsum before 
the inclusion of the Crystal Ball measurement at 3.67 GeV. The dashed curve in 
Figure 8 shows Rsum after the inclusion of the new data point. The peak of the 
between 1.00 GeV and 1.04 GeV is suppressed in both figures. 

A comparison of Figures 7 and 8 shows the effect of the Crystal Ball measure- 
ment quite clearly. Before the point is added to the analysis, there is reasonable 
agreement between the functions Raum{W) and Rfi^di^) in the region 1.0-1.8 GeV. 
Between 1.8 GeV and 3.6 GeV, Rf,^'!^^ is generally larger than Rsum- After the intro- 
duction of the Crystal Ball measurement, the 772 measurements are renormalized 
to larger values and the fitting function generally exceeds the R^^^ compilation 
throughout the region. 
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Figure 7. A comparison of our total Rhad function (dashed curve) before the 
inclusion of the Crystal Ball measurement at 3.67 GeV with that from Reference 7. 

The agreement between Rsum and -R^^"^ in the charm threshold region between 
3.6 GeV and 5.0 GeV is also quite poor. The Rsum function follows the shape of 
the DASP fit to the continuum under the ■0(4040), -0(4160), 0(4415) resonances 
and includes the resonances explicitly for comparison. The size of the continuum 
portion is determined at 3.6 GeV and 5.0 GeV by the most precise data in those 
regions (Crystal Ball data in both cases) yielding a continuous result. The R^^^ 
compilation generally exceeds Rsum throughout the region reflecting the fact that 
DASP and PLUTO generally measured large values of Rhad in th^ region. The more 
precise Crystal Ball measurements begin at 5.0 GeV and pull the R^^d function 
to smaller values, creating an apparent structure in the 4.4-5.0 GeV region. The 
apparent structure is not seen by any of the experiments that have measured the 
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shape and magnitude of Rhadi^) in this region and is entirely a consequence of 
ignoring the shape information inherent in the data (a more correct procedure 
would renormalize the data sets so that the integrated function was smooth and 
continuous in the 5-5.2 GeV region). 
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Figure 8. A comparison of our total Rhad function (dashed curve) after the 
inclusion of the Crystal Ball measurement at 3.67 GeV with that from Reference 7. 

In the region = 5 — 10 GeV, the agreement of the Rsum and R^^^ functions 
is somewhat better except for some wiggles in -R^^^ at the larger energies. Above 
b-threshold and below W — GeV (where the authors of Reference 7 begin to use 
perturbative QCD), the compilation is somewhat larger than Rsum reflecting 
the fact the the PEP/PETRA measurements of Rhad are somewhat larger than 
those predicted by perturbative QCD with currently favored values of as{M'^). 

The differences shown qualitatively in Figure 8 are quantified in Table 2 using 
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the detailed breakdown scheme presented in Reference 7. The entries in square 
brackets are from Ref. 7 before the apphcation of corrections for incomplete vacuum 
polarization correction. In Appendix A, we demonstrate that the non-application 
of this correction is generally a more accurate approximation than the use of the 
factor favored by the authors of Ref. 7. A comparison of our result with the 
bracketed quantities (or the means of the pairs of quantities) is probably the more 
relevant one. Note that our value for the contribution of 7r''"7r~ final states in the 
interval W — 0.28 — 0.81 GeV of 24.11 (in units of 10~^) is somewhat smaller than 
the value of 26.08 given in Reference 7. The difference may be due in part to the 
preference of our fit for smaller values of l-FyrP than the central values of the OLYA 
measurements between 0.6 GeV and 1.0 GeV (see Fig. 5). The opposite behavior 
is observed when the full function used in the analysis of Kinoshita, Nizic, and 
Okamoto [50] is fit to the data. The large-energy tail of this function decreases 
with energy more steeply than do the data points. A fit to this function prefers 
larger values of I-FttP than the central values of the OLYA measurements between 
0.6 GeV and 1.0 GeV yielding a contribution to Aahad(Ml) of 25.39. Excluding the 
infiuence of the steeply falling tail by restricting the fit of the KNO function to the 
region W < 1.0 GeV relaxes some of the bias and yields a Aahadi^'z) contribution 
of 24.76. These differences are reflected in the large size of the parameterization 
uncertainty given in Table 1. 

We conclude that the agreement of our analysis with one based almost en- 
tirely on trapezoidal integration is somewhat poorer than a comparison of the final 
^oihadi^'z) results would indicate. Part of the discrepancy is caused by the loss of 
shape information from multi-point measurements inherent in the averaging pro- 
cedure which treats the individual measurements as independent. An associated 
side effect is that sparse, newer measurements influence the integrated function 
only over an interval between neighboring older measurements. The addition of 
the precise Crystal Ball point (which fixes the normalization of Rhad over a large 
region in our analysis) to a trapezoidal analysis would affect only a very small 
region. Conversely, the trapezoidal analysis remains infiuenced by older measure- 
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merits until they are replaced by newer measurements at the same or very nearby 
energies. The effect of the apparent structure in the charm threshold region or 
the large Rhad values from the PEP/PETRA region will persist until replaced (or 
influenced) by newer measurements at the same energies. The use of a continuous 
fltting function in our analysis allows us to interpolate between sparse but precise 
points. For these reasons, we do indeed "believe more in the integration of our fits 
than in trapezoidal integration" as noted by the authors of Ref. 7. 

Table 2: Comparison of the various contributions to Aa^adi^'^) with 
those published in Ref. 7(in units of 10~^). The entries in square brackets 
are from Ref. 7 before the application of corrections for incomplete vacuum 
polarization correction. 



Final State 


W Interval (GeV) 


This Work 


Reference 7 


P 


0.28-0.81 


24.11 


26.08 [26.23] 




0.42-0.81 


2.87 


2.93 [2.96] 




1.00-1.04 


5.03 


5.08 [5.15] 






11.01 


11.34 [11.93] 


T 




1.18 


1.18 [1.27] 


all hadrons 


0.81-1.40 


13.55 


13.83 [13.99] 


all hadrons 


1.40-3.10 


30.42 


27.62 [28.23] 


all hadrons 


3.10-3.60 


5.62 


5.82 [5.98] 


all hadrons 


3.60-9.46 


48.16 


50.60 [50.50] 


all hadrons 


9.46-40.0 


90.67 


93.07 


all hadrons 


40.0-OO 


42.64 


42.82 


Total 




275.2 


280.4 [282.1] 
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3. Conclusions 



We have reevaluated the hadronic part of the electromagnetic vacuum expecta- 
tion value using the standard dispersion integral approach that utilizes the hadronic 
cross section measured in e+e~ experiments as input. Previous analyses arc based 
upon point-by-point trapezoidal integration which does not treat experimental er- 
rors in an optimal way. We use a technique that weights the experimental inputs by 
their stated uncertainties, includes correlations, and incorporates some refinements. 
We find the five-flavor hadronic contribution to the fractional change in the electro- 
magnetic coupling constant at = M|, Aci;(M|), to be 0.02752 ± 0.00046, which 
leads to a value of the electromagnetic couphng constant, a~^{M'^) — 128.96±0.06. 

The current generation of Z-pole asymmetry measurements have already de- 
termined the effective weak mixing angle sin^ to a precision of ±0.00028 [57]. 
Future measurements may improve the determination to the ±0.00020 level. This 
is comparable to the theoretical uncertainty of ±0.00016 which follows from the 
±0.06 uncertainty on q;~^(M|). It is clear that improved understanding of q;(M|) 
is desirable and it is also clear (from Figure 4) that improved understanding re- 
quires improved data in the 1^ = 1 — 5 GeV region. Additionally, the differences 
with the trapezoidal approach noted in Section 2.7 stem from questions dealing 
with the optimal use of rather poor quality data. Improved data will tend to make 
these issues less important. Among the active experimental programs of the world, 
only the BES Collaboration at the Beijing Electron Positron Collider is positioned 
to make improved measurements of Rhad the region W — 2 — 5 GeV. They are 
urged to include them in their long term planning. 
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4. Appendix A: Vacuum Polarization Corrections 

4.1 Corrections to Rhad 

The quantity Rhad is the ratio of s-channel cross sections and can be written 
as follows, 

where the tree-level cross sections cr°(s) are related to the physical ones (al- 
ready corrected for initial state radiation) by the simple expression, o-^(s) = 
a{s)aQ/a'^{s). Since radiative corrections calculations combine external photonic 
corrections and virtual corrections, it is more straightforward for experiments to 
extract cr^adi^) from their data than it is to extract (J}iad{s). Note that cji^(s) is 
a simple numerical constant which is applied to the measured cross section after 
radiative corrections. 

In Reference 7, Eidelman and Jcgcrlehner point out that many of the earlier 
measurements of Rhadi I-^ttP; and \Fk\'^ were corrected for leptonic vacuum po- 
larization effects but were not corrected for hadronic vacuum polarization effects. 
To rectify this problem, they make the assumption that individual experiments 
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directly measure hadronic cross sections and apply the factor, 

rf^=[l + 2Aa,(.)]^, (35) 

to all measurements of Rhad, I-^ttP) ^-iid l-^-ft'P below the and to the Mark I 
measurements below charm threshold. 

Unfortunately, the integrated luminosity for each measurement must be deter- 
mined from the measurement of an additional physical process. Thus, experiments 
rarely measure cross sections directly but nearly always measure the ratios of cross 
sections. In this case, the measured value of Rhad (or l-^P/5^/4) is determined 
from the ratio of the number of observed hadronic events Nfiad to the number of 
observed normalizing events Nnorm, 

Rhad = —J^ ^TT^T' ^ ' 

where 5rc incorporates all radiative corrections to the hadronic yield, e is the 
efficiency-acceptance product for hadronic events, and anorm is the physical cross 
section for the normalizing events (including all radiative effects) integrated over 
the acceptance used for the luminosity measurement. We note that the incomplete 
application of vacuum polarization corrections is a problem that applies to both 
the hadronic and normalizing cross sections. In this case, the actual correction 
should be 

_ aj{s) anorrni-s) ftQ + 2Aa^(s)] 

where a£(s) and crl^fyrmis) incorporate leptonic vacuum polarization corrections 
only. The difference between the two right-hand terms involves the (numerically 
insignificant) question of whether the original vacuum polarization corrections were 
performed to all orders or to first order only. 

All of the early measurements of Rhad, I-^ttP, and \Fk\^ are normahzed to the 
number of lepton pairs observed in some portion of each apparatus. Most of the 
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experiments did not have (or did not use) small-angle Bhabha scattering luminosity 
monitors but relied instead upon large-angle lepton pairs observed in the central 
region of each detector. The combination of the leptonic final states and geometric 
acceptance used by the major experiments is summarized in Table 3. Several ex- 
periments use muon pairs to normalize their results. Since the vacuum polarization 
corrections to s-channel processes can be factorized (see equation (34)), the correc- 
tion factor given by equation (37) is identically 1. The remaining experiments use 
a combination of e'^e~ and events or e^e~ events alone to normalize their 

results. The electron-pair final states are produced by the sum of s- and t-channel 
subprocesses. The vacuum polarization corrections to the dominant t-channel con- 
tributions are proportional to a^{—t). Since the t-channel contribution dominates 
the Bhabha cross section, the correction factor Tc is given roughly by the following 
expression. 



is logarithmic and the magnitude of —t at the large angles used by most of the 
experiments is comparable to s (typically, —t/s = 0.2 — > 0.4). For this —t range, 
the first ratio in equation (38) is typically a few percent less than unity and the 
second ratio is a few percent larger than unity. The net correction is therefore 
quite small. A complete calculation of the correction factor rc for requires that 
all luminosity event selection criteria be incorporated into complete calculations of 
Cnorm and CTnorm (incorporating all radiative corrections). Rather than undertake 
such an arduous procedure, we estimate the size of the correction from a simplified 
calculation which accounts for vacuum polarization effects and approximate angu- 
lar acceptance. The estimate uses the low energy parameterization of ^ahad found 
in Ref. 4. The results of this estimate are listed in Table 3 along with the correction 
advocated by the authors of Ref. 7. Note that the corrections to the pseudoscalar 
form factors are estimated assuming that the original leptonic vacuum polarization 
corrections included electron and muon contributions. The corrections to the Rhad 





a 



:2(.) ali-ty 



The key point in this discussion is that the dependence of a{q^) upon the scale 
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measurements are estimated assuming that the original corrections included only 
the electron contribution. 



Table 3: Summary of the incomplete vacuum polarization correction fac- 
tor Tc and that of Ref. 7, r^*^. 



T ^ , , 

Hixp. 


Meas. 


Norm. 




cos(^)| 


W (vjeVj 






NA7 [13] 






— .8(0— ^.yy/ ' 


0.320 


1.0000 


0.9982 










0.422 


1.0000 


0.9972 


OLYA [12,20] 


\F.\\ \Fk? 


ee + yU/x 


<U. / i 


0.400 


0.9984 


0.9974 










1.397 


0.9952 


0.9893 


CMD [12,21] 


\FA\ \Fk? 


ee + //// 


<U.DU 


0.360 


0.9988 


0.9978 










0.820 


0.9970 


0.9934 


TOF [14] 


\F |2 


ee + //// 


<U.Z4 


0.400 


0.9990 


0.9974 










0.460 


0.9988 


0.9968 


^X7: [17] 


IF |2 


ee 


<0.61 


1.250 


0.9958 


0.9902 










1.520 


0.9955 


0.9886 


MEA [19] 


77' Z 77' Z 

FttI , \i'K\ 


ee 


<0.77 


1.6 


0.9941 


0.9826 






fjjfjj 




1 

1 .^o 


l.UUUU 


u.yooo 


DM1 [16,22] 


WA^. \Fk? 


ee 


<0.50 


0.480 


0.9983 


0.9966 










2.060 


0.9960 


0.9860 


DM2 [18,23] 


\FA\ \Fk? 




<0.87 


1.350 


1.0000 


0.9896 










2.400 


1.0000 


0.9848 


772 [26] 


Rhad 


ee 


<0.64 


1.42 


0.9933 


0.9839 










3.09 


0.9935 


0.9757 


Mark I [28] 


Rhad 


ee 


<0.60 


2.60 


0.9936 


0.9772 










3.65 


0.9958 


0.9756 


DASP [29] 


Rhad 


ee 


<0.7l('^) 


3.6 


0.9946 


1.0000 


PLUTO [30] 


Rhad 


ee 


0.98 


16^0.9977 


3.6 


0.9756 


1.0000 


CMD [53] 




ee + /i/x 


<0.60 


0.782 


0.9971 


0.9904 


ND [54] 




ee 


<0.65 


0.782 


0.9942 


0.9904 



Interval in cos^. 

*^^-*Used small-angle e+e~ events normalized to this large angle region. 



The reader should note several things. The corrections to the Rhad^ I-^ttP, and 
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\Fk\ measurements are always a factor of seven or more smaller than system- 
atic normalization uncertainties associated with the measurements. In all cases, 
the correction applied by the authors of Ref. 7 overestimates the true size of the 
correction. This overestimate is small where the correction is small but becomes 
significant at larger energies where the Eidelman-Jegerlehner correction exceeds 
1%. In this region, the non- application of the correction (rc = 1.0) is a better ap- 
proximation than the one used by the authors of Ref. 7. The Eidelman-Jegerlehner 
analysis did not correct the hadronic continuum measurements of the DASP and 
PLUTO Collaborations at charm threshold although it appears that neither group 
applied hadronic vacuum polarization corrections [58]. The normalization DASP 
measurements was determined from the total number of large-angle Bhabha scat- 
tering events and is subject to a small correction. The PLUTO experiment nor- 
malized its measurements with a small-angle luminosity monitor which sampled 
a region of small —t. The cancellation of the vacuum polarization corrections is 
correspondingly smaller and the correction is larger. 

4.2 Corrections to Resonance Parameters 

The Breit-Wigner cross section used in Section 2.5 to calculate the resonant 
contribution to Aa^adiM'^) requires the mass, total width, and electronic width 
of each resonance as input. The electronic widths Tee are defined to be physical 
quantities (not corrected for vacuum polarization effects) and differ from the tree- 
level quantities Fgg that have been used often in the past. The electronic widths 
for narrow and broad resonances are determined by different techniques but are 
always proportional to the peak hadronic cross section of the resonance (measured 
in e'^e~ collisions) or to the measured energy-integral of the hadronic cross section 
(taken over the resonance), 

where all quantities are defined in equation (36) except for 5'^.^. which accounts 
for radiative corrections to the hadronic yield but excludes vacuum polarization 
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corrections. The inclusion of vaccum polarization corrections into (5^^ — >■ 5rc) 
yields a measurement of the tree- level quantity Fgg. 

As in the case of the cross section and form factor measurements, many of the 
older measurements of the electronic widths were not corrected for hadronic vac- 
uum polarization effects. It is clear that measurements of Fgg must be corrected by 
the same correction factor rc defined in equation (37). However, for measurements 
of Fee, vacuum polarization corrections to the hadronic yield are not applied and 
the appropriate correction factor Qc pertains to the normalizing cross section only. 



^normi^) 
Is)' 



9c = (40) 



As was discussed in Section 2.5, the Review of Particle Properties hsts phys- 
ical widths for the ip- and T-family resonances as derived from their own fitting 
procedure. The electronic width of the 0(1020) is either the physical value or an 
average of the tree-level and physical values and is assumed to the the physical one. 
The oldest measurements of these quantities were corrected for electron vacuum 
polarization effects only and require the application of the additional correction 
factor Qc- Estimates of this factor are listed in Table 4 for measurements of the 0, 
J/iJj{lS), and J/ip{2S) electronic widths. The weighted average of the (p correc- 
tion factors is applied to the PDG value of Fee(0)- The corrections to the z/'-family 
measurements are quite small if the original measurement was normalized to small- 
angle Bhabha scattering and can be as large as 2% if the large angle cross section 
was used as a normalization. Unfortunately, since the quoted electronic widths are 
derived from global fits, it is difficult to estimate the effect on the final value of Fee- 
Therefore, we do not apply any corrections to the electronic widths of the ■0-family 
but we do inflate the uncertainties on Fee by the size of the largest correction. 

Unlike the other resonances, the electronic width of the a; (782) listed in the 
Review of Particle Properties is the tree- level one. We therefore apply the weighted 
average of the correction factors rc hsted in Table 3 for the dominant CMD and 
ND measurements. 
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Table 4: Summary of the incomplete vacuum polarization correction fac- 
tor Qc. 



Exp. 


Res. 


Norm. 


cos(^)| 


9c 


DM1 [55] 


0(1020) 


ee 


<0.50 


1.0071 


OLYA [56] 


(^(1020) 


ee 


<0.71 


1.0052 


Mark I [59] 


jms) 


ee 


0.9997^0.9999 


1.0000 


772 [60] 


jms) 


ee 


0.9945^0.9986 


1.0002 


MEA [61] 


jms) 


ee 


<0.77 


1.0158 


DASP [62] 


jms) 


ee 


<0.71 


1.0169 


Mark I [63] 


J/^{2S) 


ee 


<0.60 


1.0204 


DASP [62] 




ee 


<0.71 


1.0189 
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